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Abstract 

Motivated by the Kaluza-Klein theory with large extra spacetime dimensions the neutrino star 
built from the massive sterile neutrinos core and the massless brane neutrinos envelope is presented. 
The six-dimensional compactification scale ~ 15 MeV gives maximal neutrino mass M max = 2.3 x 
lO 4 M with radius R 

max — 1-2 x 10^ km. The maximal neutrino star parameters varies with 
temperature. In the limit of the neutrino ball approximation the maximal sterile neutrino star is 
M max = 1.1 x 1O 6 M . 
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Introduction 



Recently there has been considerable interest in the field theories with large extra space- 
time dimensions. In the comparison to the standard Kaluza-Klein theory these extra di- 
mensions may be restricted only to the gravity sector of theory while the Standard Model 
(SM) fields are assumed to be localized on the 4-dimensional spacetime This is 

a promising scenario from the phenomenological point of view because it shift the energy 
scale of unification from 10 19 GeV to 1 — 100 TeV. It has been recently shown Q that 
this framework can be embedded into string models, where the fundamental Planck scale 
can be identified with the string scale which could be as low as the weak scale. The extra 
dimensions have the potential to lower the unification scale as well ||. 

The aim of this paper is to examine the degenerate neutrino star originated from the 
extra dimensional theory. The neutrino star ( neutrino ball ) was a subject of interest in 
the theory of the Standard Model || . The fermion star in the extra dimensional theory was 
also the a subject of interest ||.The neutrino star model is capable to explain the nature of 
the object in Sgr A* in the center of the Galaxy 0. 



The bulk neutrino extension of the electroweak theory 

Much of the interesting phenomenology of brane-world models is associated with the 
Kaluza-Klein theories || that originate from large, gravity-only additional dimensions. The 
higher-dimensional bulk fermions lead to Kaluza-Klein towers of standard model singlets that 
may be interpreted as sterile neutrinos [JR], O, [T2], [IB]. In this paper we shall consider the 



six-dimensional Kaluza - Klein theory. Let us now consider the action in the six- dimensional 
spacetime: 

S = J d 6 xV=f 6 L = S g + S F = J d 6 x^/=g~ 6 {L g + L F ), (1) 

where g% = det(gMN) and M = {/i, i}, N = {v, j} with x M = {x^.y 1 }, i — 1,2. The 
metrical tensor in the six-dimensional spacetime can be written: 

I e- 2 ^g, u 
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y -^. e + 2 ^)//o 



(2) 



According to the above definition we can write: 

V^=V^e~^». (3) 



We consider here the Lagrangian of the field as follows: 

L = L g + L F , (4) 

Lp = L brane 5{y) + L bu i k , (6) 

where kq is the six-dimensional gravitational coupling. Its natural interpretation originates 
from the distance scaling in the four-dimensional spacetime. Let us compactify the six- 
dimensional spacetime (M.§ — > M.^ x5'x S 1 ) to the four-dimensional Minkowski one. In 
this paper we assume that the extra dimensions are compactified on the 2 dimensional torus 
with a single radius r 2 . The six-dimensional action may be rewritten as: 

S = J d A x J d 2 y^—g e L = J d 4 xy/—gC, (7) 

where / d 2 y = (27rr 2 ) 2 . The six-dimensional gravitional coupling = 87rG6 is convenient 
to define as 

where M - is the energy scale of the compactification (~ 10 — 100 TeV). Cosmological 



consideration |T4| gives the bound M > 100 TeV what corresponds r 2 < 5.1 x 10 mm. If 
we define the four- dimensional coupling constant k = 8tiGn = 8nMp 2 we get 

M Pl = ATiM A r 2 2 . (8) 

The parameter f is defined as: 

ft = -M 2 Pl (9) 

7T 

to produce the 1/2 term in for the dilaton field kinetic term. The model can be motivated 
by the ten-dimensional string theory with the string scale M s . For simplicity we assume that 
the first 2 dimensions have the same size characterized by the radius r 2 . This is likewise 
assumed for the remaining 4 dimensions where the corresponding radius is called < r 2 . 
The four-dimensional Planck constant 

( 1F )2 = 167r ( 2 ^M s ) 2 (27rr 4 M s ) 4 (10) 
is connected to the string scale M s . In the result of the comparison to (§) we have 

M = v / 47r(27rr 4 M s )M s . (11) 



M s 


M 2 = r^ 1 


r-, 1 


M 


1.2 TeV 


1.6 MeV 


432.05 MeV 


7.42 10 7 GeV 


1.2 TeV 


15 MeV 


141.11 MeV 


2.2710 s GeV 



TABLE I: The paremeters set of the model [15] 



The paremeters of the model |15[ are presented on the Table |. In this section we shall 
extend the Standard Model minimally with the bulk neutrino ip(x,y). The lagrangian of 
the neutrino sector of the model is then: 

Lbrane = iLfT^D^Lf - — (/l/ )0 ^ a L/#Q + h.c), 



and the fermion field are 



/,.o 



V 



f 





/ = e, v, t. 



(12) 



The Higgs field will have only the residual form 

H = 




(13) 



after the spontaneous symmetry breaking. The additional bulk neutrino is described by the 
Lagrange function 



£buik = #r d N ip -m D (ip L ip R + h.c), 



(14) 



where ipn t L = (I ±T 7 )ip. The system has global U(l) lepton symmetry which generates the 
lepton charge Ql- The bulk neutrino field can be decomposed into four-dimensional Dirac 
spinors t/> , where a = {+,—}• The gamma matrices T are defined as 



MM i 



{r M ,r«} = 2 9 

In the flat four-dimensional spacetime, when g = r]^ u they may be defined as 



'(2), 



(15) 



(16) 
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and 



{y\y} = 2rr/, 
{ 7 l ,7 J } = -2^ 
4 



where 



r 7 = r°...r 5 = 7 5 ®a 3 , 7 5 = -z 7 ° 7 V7 3 - 



Using the metric tensor six-dimensional form (y) and the ( |16|) we can calculate the Lagrange 
function in the four- dimensional Minkowski spacetime 

L = (LfD^L + e- 2f//o J2 t^aJx)YD^ a ,n(x) (17) 

-m D 1Pa,n,L^a,n,R ~ e €// ° ^ ^a,M (Y<)^a,n (x) , 

where 



a. n 



< = — n *> ml=yl = m^{n\ + nl). (18) 



1 i 2 2 1 

-n, m„=^ = m- 
The bulk neutrino is decomposed as 



e ?//o ]_ 

i } a{x,y) = 7 - zY^a,n(x)exp(—iny). (19) 



Each of these four-dimensional Dirac spinors can be decomposed into Weyl spinors 

/ 



(20) 



\ Va,i 

After electroweak symmetry breaking we introduce the mass 

mfa = h /' aV „ 10- 1 MeV. (21) 
f ' a 27r(M s r 2 ) 1 ; 

The exact diagonalization of the mass matrix gives three exactly massles Weyl fermions 

= U f,f'4 + Vf,a(0)Va(0) (22) 

and two Dirac spinors for each mode number n with mass 

Mn = sfml + < (23) 



In general ||15| , there is a superposition of the electroweak neutrinos on the brane and 



the Kaluza-Klein bulk neutrinos t] a (n). 



The neutrino star 



An alternative model for the supermassive compact object in the center of our Galaxy 
has been recently proposed by Tsiklauri and Viollier [plj[|17|]. The main ingredient of the 
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proposal is that the dark matter at the center of the galaxy is non-baryonic composed with 
massive neutrinos or gravitinos. Such neutrino balls could have formed in early epochs, 
during a first-order phase transition in the Standard Model. 

In case of the spherically symmertic gravitational field we have: 

\ 

_ e X(x) 

,2 

— r 2 sin 2 6 



( e "W 
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\ 



(24) 



/ 



In the similar way the dilaton filed £(r) will be dependent on the radius r. As £(r) ~ 
fo = MpxjypK ~ 10 19 GeV A we shall neglect the dilaton filed £(r) in the first approximation. 
In this paper we present numerical results describing the structure of neutrino star. It is 
possible to describe a static spherical star solving the Oppenheimer-Tolman-Volkoff (OTV) 
equations (more general case with the dilaton filed £(r) is presented in Appendix I). 



dP(r) 
dr 



G 



(p(r) 



dm(r) 
dr 



P(r)Am(r) + ^P(r)r 



2Gm(r) 
c 2 r 



<iirr 2 p(r). 



(25) 



(26) 



Having solved the OTV equation the pressure p(r), mass m(r) and density p(r) were ob- 
tained. To obtain the total radius R of the star the fulfillment of the condition p(R) = 
is necessary. Our aim is to achieve the equation of state of neutrino star matter at finite 
temperature. In such a case the physical system can be defined by the thermodynamic 



potential |18 



n 



-k B TlnTr(e-^ H -^ QL) ), (27) 

where (3 = 1/ksT, is the Boltzmann constant, Ql lepton charge. The chemical potential 
p reflects the lepton number conservation. H stands for the Hamiltonian of the physical 
system. All needed averages are calculated with the Hamiltonian H. We define the density 
of energy and pressure by the energy - momentum tensor 



T^v = (P + t)u^u u - Pg 



where is a unite vector = 1). So, the calculations give 

e(xp, T) = pC 2 = epfiulk + ^F,brane, 



(28) 



(29) 



P(x F ,T) = P F , 

bulk ' ± F branei 

(30) 

where 

£f = e x(x F ,T), (31) 

P F = P o 0(x F ,T). (32) 

The fact that the bulk Dirac neutrinos are massive means that they play the same role like 
ions in a white dwarf or neutrons in a neutron star. We have 



"X.bulk 



| poo | 

(x F ,T) = -s dzz 2 V^Tl{ 7 r 7T =^ /w r— (33) 

77 Jo exn( (\/l 4- — ii ) t) 4- 1 



+ 



•exp((v / T + ^-/i , )/r) + 1 
1 

exp(( v / rTl2 + /i')/r) 4-1' 

1 /*°° r 1 

W(x^T) = — y Q ^_ { exp( ^ i + ^_^ ) + i (34) 

1 

+ exp((v / rT^ + / u')/r) 4-1 



where r = (k B T)/M< 



2- 



fj/ = fj,/M 2 = VI 4- 4 (35) 

and 

x F = A; F /M S . (36) 



Similarly to the paper JT9[ we have introduced (p5| , p6D the dimensionless "Fermi" momen- 
tum even at finite temperature which exactly corresponds to the Fermi momentum at zero 
temperature. For the massless brane neutrinos we have 



Xbrane(x F , T) - ^ 2 J dz Z 3 { ^ {z _ Xf)/t + ^ + ^ {z+Xf )/ t + ^ } ( 3 ?) 

r eiz -tyl + 1 {^(^ /T ) + ^(^ /T )} , (38) 



where the presure is 



^>brane(x F , T) — -tbrane{.X p , T) . (39) 



Both e F and Pp depend on the neutrino chemical potential fi or Fermi momentum x F . This 
parametric dependence on \i (or xp) defines the equation of state. 
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Similarly to the paper [HJ we have introduced the dimensionless 'Fermi' momentum even 
at finite temperature which exactly corresponds to the Fermi momentum at zero tempera- 
ture. Both ep and Pp depend on the neutrino chemical potential /i or Fermi momentum xp. 
This parametric dependence on fi (or xp) defines the equation of state. 

When Fermi momentum reaches the second Kaluza-Klein level ~ the number of 

available neutrino modes will change. For the high bulk neutrino density in the ultrarela- 
tivistic limit all Kaluza-Klein modes should be included. In that limit energy density of the 
bulk neutrinos is equal to 



Xbuik(xF, T) = ^(r 2 M 2 )Jdzz*{ {dz _J )/T + 1) + {e(z+XF ) /T + 1) } (40) 
(r 2 M 2 )t«f Q e{z _ XF)/T + i {U 6 (e^) + ^(e-^)}- (41) 



7T 

In the ultrarelativistic limit the equation of state is 

4>buik(xF,T) = -e bu i k (x F ,T). (42) 

The equations ( p5| , p6|) are easy integrated numerically, For example, for the neutron star 
with the central density p c = 10 4 g/cm 3 and temperature T = 50keV the star density and 
pressure profile are presented on the Fig.[TJ and Fig. 0. Similarly to a structure of white 
dwarf massive sterile neutrinos like ions contribute to the density of the star , while the 
massles neutrinos like electrons contribute to the pressure. This feature will be more visible 
with the increasing temperature. The parameters of the maximum mass configuration are: 

M max = 2.3 x 1O 4 M , R = 1.2 x 10 6 km. (43) 

This fact is easy to notice on the mass-radius diagram (Fig. ||D . In general, we have all family 
of neutrino stars, depending on growing neutrino Fermi momentum. This result concerns the 
Ferni momentum below the second Kaluza-Klein level ~ mp,\[2 . The ulrarelativistic limit 
gives higher masses M max ~ 1.6 x 1O 6 M . The density profile in this limit is presented on 
the Fig[l| (magenta). In limit of bulk neutrino ball (Appendix I) the neutrino star properties 
are presented in Table [TT]. and the neutrino star mass dependence from the central density 
Pc (Fig. §. 
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FIG. 1: The density profile of the neutrino star for the bulk neutrino (blue), brane ones (red), in 
the ultrarelativistic limit (magenta) and in the neutrino ball case (the dot line). 




FIG. 2: The pressure profile of the neutrino star. 

Appendix 

The Einstein equations in in the six-dimensional spacetime (the metric tensor (@,0) ) 
can be written as 

+ 1 _ 2 PM!) _ L = ^ e -™.p W , (44) 

r r l Jq r z c 4 



e -»^.m + L + a^LL) _ I = JJIp e -^m p{r) , (45 ) 

r r z Jq r l c z 
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FIG. 3: The R-M diagram for the neutrino star. 
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FIG. 4: The neutrino star mass dependence from the central density 
The Einstein equations give also the equation for the dilaton field: 

C'(r) + fcf (r) - ±f (r)(A'(r) - v'{r)) - ^' (rf+ 
HU^ r) - !) + U X> (r) - «/(r)) + §(A'(r) - i/(r))i/(r)- 
XM + ^e A M- 2 « r » P s (r)). 

The continuity equation T.^ N = gives 

1 1 (P + C*p) f (P + CV) ' 
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We assume that the energy-momentum tensor has a diagonal form diag{T^ v ) = 
{e, P, P, P, P s , P s }. Integrating of equation ( flop yields we can write: 

p -x(r) = /^((l-^p^e fo )e J "4, ) 

exp(/dr(± + |£'(r) 2 )) " 1 ) 



If we define 



e -A(r) = x _ f^» m ( r ) (49) 



then we can obtain the generalized Oppenheimer-Tolman-Volkoff equation: 
dP(r) G 6 (p(r) + 3£)(m(r) + f P(r)r 3 e ~ 2 «M//o) 



J/- /-2 a -f£m(D) f (50) 

c 2 p + P). (51) 



X'{rf £'(r)(P-2P s -c 2 p) w 2 



v / 2 A (P + c 2 p) 
In vacuum P = , P s = 0, of course we have well known Schwarzschild solution 



e -A(r) = (1 _ 

r 

g „(r) = (1 _ !£) 

r 



£(r) = 0. 

Neglecting the dilaton field £(r) the Oppenheimer-Tolman-Volkoff (^5[) . 

The theory of neutrinos, bound by gravity, can be easily sketched considering a Thomas- 



Fermi model for fermions [p0| . We can set the Fermi energy equal to the gravitational 
potential which binds the system, and see that the number density is a function of the 
gravitational potential. Such a gravitational potential will obey a Poisson equation, where 
neutrinos (and anti-neutrinos) are the source term. Including gravity the local equilibrium 
condition demands 

7^W)Tr ] = - 4lG " p(r) ' (K) 

Inside the ball the pressure of the massive bulk neutrinos is P = p/5 while P = p/3 for the 
brane neutrinos in the relativistic limit in high temperature. Defining 

p(r) = p c e* (53) 

with 

« = 24nG N p c (54) 



11 



Pc (g/cm 3 ) 


R (km) 


Af(M ) 


10 3 


5.9 10 6 


1.1 10 6 


10 4 


1.8 10 6 


3.5 10 5 


10 5 


5.9 10 5 


1.1 10 5 



TABLE II: The bulk neutrino star properties in the neutrino ball approximation, 
we have the Liuoville equation 



The Laplace operator is 



1 d , 2 d , 
r 2 dr dr 



Using the thin-wall approximation one can obtain the following expression 

F (r) = ' 



cosh 2 (4^)' 



(55) 
(56) 

(57) 



The equation (57) allows to estimate the mass of the neutrino ball which is given by the 
dependence 

roc dw 3 

M(R) = 4vrp c / = (58) 

with the radius scale 

= J%. (59) 



The bulk neutrino density profile (|57D is presented on the Fig. |l] (the dotted line). 
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